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Abstract 


We deal with some effective numerical methods for solving a class of 
nonlinear singular two-point boundary value Fredholm integro-differential 
equations. Using an appropriate interpolation and a q-order quadrature 
rule of integration, the original problem will be approximated by the non- 
linear finite difference equations and so reduced to a nonlinear algebraic 
system that can be simply implemented. The convergence properties of the 
proposed method are discussed, and it is proved that its convergence order 


will be of O(nmint.9-3h), Ample numerical results are addressed to con- 
firm the expected convergence order as well as the accuracy and efficiency 
of the proposed method. 
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1 Introduction 


In this study, we consider the following nonlinear singular two-point boundary 
value Fredholm integro-differential equation (SFIDE): 


1 


(t%y'(t))’ = f(t) +f u(t, s)u(y(s))ds, t€(0,1], O<a<1, (1) 


y(0) =a, yl) =6, (2) 
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where f(t), y(t) and kernels v(s, ¢), u are known Lz functions, and all of them 
are in C*4((0,1]). The nonlinear singular problems are extensively arisen in 
many applications in physics and astrophysics [11, 9, 12, 19, 20], chemical 
and mechanical engineering [18, 5, 15, 13], physiological process [14], popula- 
tion dynamics and epidemiology [23], fluid mechanics, electro hydrodynamics, 
nuclear physics, and chemical kinetics [2, 10]. 

Also, the nonlinear singular problems have many applications in stel- 
lar structure, thermal explosions, isothermal gas spheres, radiative cooling, 
thermionic currents, and the thermal behavior of a spherical cloud of gas 
[3, 17, 22]. 

The majority of engineering applications and various branches of sci- 
ence, such as financial mathematics, oceanography, population dynamics, 
fluid mechanics, plasma physics, electromagnetic theory, artificial neural net- 
works, and biological processes, have been dominated by Fredholm integro- 
differential equations (FIDEs) [4, 6]. In general, the analytical solution of 
FIDEs is not available. As a result, various numerical techniques for deter- 
mining approximate solutions of FIDEs have been created. The situation 
is significantly more complicated for FIDEs with singularities. A particular 
type of them called singularly perturbed Fredholm integro-differential equa- 
tions (SPFIDEs), was discussed in [1, 6, 8, 7]. Numerical analysis of SPFIDEs 
has not yet been widely utilized. In this study, we focus on a specific case 
of nonlinear singular two-point boundary value Fredholm integro-differential 
equations of the form (1)-(2). Since solving problems of this type is very 
difficult, the main motivation of this study is to construct an efficient and 
useful numerical method with O(h2) accuracy in the Ly norm for nonlinear 
singular problems of the form (1). 

To formulate some accurate and effective methods for (1), we first ap- 
ply a finite difference method to discretize the singular ordinary differential 
equation part and a suitable quadrature rule of integration for the singular 
two-point boundary value Fredholm integro-differential part of (1). 

Then, the original problem is converted into a system of nonlinear alge- 
braic equations. The numerical solution of the derived nonlinear system is 
computed by using some solver like the Newton method. Also, the conver- 
gence analysis of the present method is established. 

The main features of the new method are as follows: 


e It can be simply implemented by converting the singular problem into 
a system of nonlinear algebraic equations. 


¢ The convergence rate of the proposed method is O(h*) with respect to 
the L.. norm when applied to nonlinear singular problems. 


e The proposed method is successful in solving some classes of singular 
problems, such as SFIDEs and SPFIDE. 
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e The provided comparative numerical simulations confirm that the pro- 
posed method is more accurate than the existing methods reported in 
the literature. 


2 Formulation of the method 


In this section, we formulate a novel numerical method for solving the two- 
point boundary value Fredholm integro-differential equation (1). At first if 
one takes q(t) = t*y’(t), then (1) reduces to 


J) =F)+ [ v(t,s)u(u(s))as (3) 


Consider the partition {t, = kh : k =0,1,...,N} ofthe interval [0,1], where 
to = Oandty = landh= x denotes the step size. For k = 0,1,..., N, let Y; 
and Vz, denote the approximate values of yx, := y(t,) and Vin := U(tK, tn), 
respectively. For (3), we can conclude that 


a)—a= fH jag ff oe,syuty (s))dsdé. (4) 


Dividing both sides of (4) by t* and then integrating over [t,,t,41] and 
[th— itr], we have 


i (v@-) a= f aca f [{ W(E5) a y(s))dsdédt 


By changing the order of integration, we get 


il 
ws — we FUT iags = fe + f wluls))uE (eas (5) 
where 
we-dcg psn gees 
Tr-1 = ; d 
pie bi 4 i 1 = #(6)dé, 


and k=1,...,N —1. Eliminating q, in (5) concludes that 


7 1 
(Ye —Yr-1 t+ f, ) + Tr (Ye — Yeri t+ fie) 
1 is 


+f utu(ey (atv (0) + Lug (sy) de =0.(7) 
[ Ca 1 Th ) 


Iran. j. numer. anal. optim., Vol. 13, No. 3, 2023,pp 444-459 


4A7 Effective numerical methods for nonlinear singular... 


To solve (7), it is sufficient to utilize some suitable numerical integration 
methods to approximate f+, v*, and its integral part. By using the in- 
terpolating polynomials of f(&) and v(€,-) at nodes t, and trai, we can 
approximate the integrals given in (6) as follows: 


fe = G0 ~F (te) + OTF ther) + 02. (the) + a3 FEE) 
aie ele 2 al 3 

Uy (s) = ag ,U(th, 8) + ay ,v(teti, s) an as, Se (t, 8) tt a 5 5, See (t, 8) t=c+’ 
Str =S;, 


in which £7 ,¢, € (te—1, tx) and &, G+ € (te, te41) and 
gy Daeg pee ygs 
or t —a-j  42-a— 

Gok = Fwy (i) kt1 th, ti 


0 
sg up (uri -o gy 


2h =e a-j\W 
,=+5 » po (waren, ty 
Che Gey) etat -m 
- 3 (ar - yg, 8 
= ay po (Huard. ay 


Assume that the functions f(t) and ot s) are continuous. Then there 


are the values cx, © (te—1, te41) and Cx, a € (min(€;), max(&;)) such that 


1 1 


fe tafe = ORF (te) + 4% f(te—1) + oft F(te1) + ex (Sf); (9a) 
Tiel Tr 
uz (s) + uf (5) = p2v(te, 8) + pp vlte—1,5) + vt Utes1, 8) + ex(v(-,5)), (9b) 
Leni Tr 


(a 2 Ste 1 ae 1 
where Pre= bo,k = 7202,k5 Pr = Te 1,k + 72 02,k: 


en(f) = —Fph?bo,nf (ce) + d3,nf (Ce), 


4 3 
ex(v(-,8)) = hee E(t s)| + beet )| 
t=Cr t=Cr 


and bj, = HG + aot 1 = 0,2,3. Finally, by applying the relations (9) 
and utilizing a suitable numerical quadrature method of order q with weights 


w = (wo, w1,-..,wy)', (7) can be reformulated as 
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N 
1 1 1 1 . ‘ 
14 ; +P, th) Vin =e + O(hit), 
ida ( x) 7, vet + Fe Zee kn = €r(f) ( ) 


(10) 
where fork = 1,...,N —landn=0,1,...,N, we set 
Un = U(Yn), 


F, = p? fre + Py fk—1 + Be TKaAs 
Vin = PrUE.n + Dy VE-1,n + Ce eee 


Therefore, an approximate method to solve the problem (1) can be formulated 
as follows: 


1 1 1 1 aul ‘ 
ee Y, Yeti ta we S0. 
( Ge x) a ae a ae o 


(11) 
where U,, = u(Y;,). Take note that the Newton method can be used to solve 
the derived nonlinear equations. Let us set Y = (Yj,...,Yn—-1)', F = 


: . P : i T 
(Fi, seey Fy_1)',V = lv. waits WN) ; and ye = (Vis Va.k; seey Vn-1,k 
fork =1,...,N—1. Then the matrix formulation of the proposed method 
(11) is also written in the following form: 


TY+AL u(¥) = —F— woloV° — wnUnV™% + 7¥oli + twY¥nIn_1, (12) 


where W = diag(wi,...,wn_—1), L = VW and 


T = tridiag(—Ty, eT o+Ty »—-Ty 2); (13) 


is a tridiagonal matrix with TX = [7,,7,41,.-- Trl) and tT, = 7: for k = 
0,1,...,N. The symbol I; signifies an (NV — 1)-column vector with entry 1 
in the position 7 and 0 elsewhere, where i = 1, N — 1. 


Remark 1. It is worth noting that according to the relation t¢y” = (t®y')'— 
at*—1y’, the presented technique may be utilized for the singularly perturbed 
Fredholm integro-differential equations discussed in [1, 6, 8, 7], as well as the 
singularly perturbed boundary value problems considered in [10, 13]. 


2.1 Convergence analysis 


In this section, the convergence analysis of the presented method (11) to solve 
the SFIDE (1) is performed. To this end, we set y = (y1,y2,---;yn—1)'- 
Then a matrix formulation of (10) is also derived as 
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Ty+hAL uly) =E(f) — F— woUoV° — wnUnV™ 4+ 0Yoli 
+ tvYnIny-1 + O(h*!)1y-1, (14) 
where E(f) = (e1(f), e2(f),---,en-1(f))' and 1y_-1 = (1,1,...,1)' € 
R‘-1. By subtracting (12) from (14), the error equation can be derived 


as - 
T(y — Y) +h L(u(y) — u(Y)) =E, 


where E = E(f) + O(h?+!)1n_1. Thus we have 
(T+hLJy)(y — Y) =E, (15) 


where Jy is a diagonal matrix containing the Jacobian of kernel u(y). That 
is, Ju = diag([Zu(y)|,_, Jp )- To formulate an upper bound for the 


L2-error ||y — Y4|| derived in (15), we first prove the following lemmas. 


Lemma 1. Let A € R”"*” and x € R”. Then 
|| Ax|]2 2 Omin(A)||Xll2, 
in which oyjn(A) is the smallest singular value of matrix A. 


Proof. If we consider the singular value decomposition of the form A = 
SZ", where S and Z are orthogonal and © = diag(o1,...,0n) is the diago- 
nal matrix with singular values o,,k = 1,...,n. Then from the orthogonality 
of S and Z, we have 


Axl, = |SE2"x| = [SZ x)],, = [22"x 
O71 an AA ii n n 2 
~ : =| a (>: su] 
On ar en ivi 2 hot =r 
where x = (#1,...,2y)'. Let omin(A) = min{ox;k = 1,...,n}. This con- 


cludes that 


n n 2 
|| Ax||, > Omin(A) ( su = Omin(A)||Z"x||2 = Omin(A)||x||2- 
k 1 


=] {= 


In the following, we may try to construct lower and upper triangular 
matrices L; and U, such that the tridiagonal matrix A can be expressed as 
the product A = L,U, of the form 
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dy Cy ] 
Cc 
be dy C2 EB ] Plc 
fy 1 p2 
An—2 dn—1 Cn-1 “ie Es "a Cn-1 
aAn-1 dn In—1 1 Pn 
a, ee se No, 


————S 
A Li an (16) 


Indeed, multiplying DL; by Uj, yields the following recursive relations: 


a 
pi = dh, k= Proi = dr41—lece, &=1,2,...,n—1. (17) 


In the following lemma, we will exhibit that the decomposition (16)—(17) for 
the matrices with the property of strictly diagonally dominant is unique. 


Lemma 2. If A is a strictly diagonally dominant matrix, then it has a unique 
LU-factorization of the form (16)—(17). 


Proof. It is sufficient to show that the elements pz, in (17) are nonzero for 
k =1,...,n. It can be done by induction. So, we show that |pz| > dx + |cel, 
where 0, = |dg| — |ag—1| — |cx| > 0,& = 1,2,...,n, and ap = cp, = 0. Since 
\pi| = |di| = 61 + |ci|, assuming |p,| > 6, + |cx| for some k = 1,2,...,n—1, 
concludes that aa < 1. Therefore, according to (17) and the strictly diagonal 
dominant of the matrix A, we have 


|ax||cx| 


QAkCk 
z \di+1| = pee| 


\Pe+il = |de+i — lece| = |degi — 


> |de+4il — lax| = beta + len+il- 


In the next lemma, we can see that under some conditions there is an 
LU-factorization in the form (16)—-(17) for every weak dominant tridiagonal 
matrix. 


Lemma 3. Assume that A is a tridiagonal matrix with the property of 
weakly diagonally dominant. If in addition |di| > |ci| and a, 4 0,k = 
1,2,...,n —2, then it has a unique LU-factorization of the form (16)—(17). 
Moreover, if d,, #0, then A is nonsingular. 


Proof. The proof of this lemma is similar to Lemma 2. 


In the next step, we investigate some properties of the matrix T given by 
(13). Since T is a weak diagonal dominant symmetric matrix with positive 
diagonal elements, then it is a positive semidefinite matrix. It is easily seen 
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that the elements of T° and T} are not vanished. Consequently, from Lemma 
3, we conclude that |T| 4 0 and p, > 0 for k = 1,2,..., N—1. This yields that 
all of the singular values oj as well as eigenvalues A; of T are positive. Now we 
establish that Amin(T) > h. To this end, we first define M = iT —I, in which 
I is the identity matrix of order N — 1. So, if A is an eigenvalue of T, then 
(FA — 1) is the eigenvalue of M. Hence, it is sufficient to prove that Amin(M) 
(the smallest eigenvalue of M) is positive. Since M is a strictly diagonally 
dominant matrix, then from Lemma 2, there exists a unique LU-factorization 
of the form (16)—(17) for this matrix with the following coefficients: 


ees a) = 1, k=1,2,..., 


Remark 2. Since h = %, the coefficients given by (8)-(9) can be approxi- 
mated as 


to 2 3 5 


h 
k ote -—a 
Ah? ak 6 Fe? 0,k , 2,k 12’ 3,k 54° 


(18) 


as h + 0. Therefore, the elements y;, yf, and V” of the matrix V will be of 
order O(h). Finally, we conclude that the elements of the matrix L will be 
of order O(h). 


Lemma 4. If there is an LU-factorization for the matrix M in the form 
(16)—(17), then |M| > 0 and p, > 0 for k=1,...,N—-1. 


Proof. Putting py = hpr, I, = —ly and utilizing (17) yield pry1 = Tesi t+Te— 
h — ply. Therefore, we get, 


Ve+1 
Yes t1—i, — Af te’ 


[koi = 


in which 7441 = ,k = 1,...,N—1. It should be mentioned that 1; = 


TT ? 
aan < 1 and im, ik = 1. So, for sufficiently small h, assuming J, < 


= we have0 <i <1,k= 
Types oN — 1. Therefore, with the help of (18), we get 


Prti > Troi -h>k. 


From Lemma 4, we can reach that the determinant of all upper-left sub- 
matrices of M is positive. Thanks [24, Theorem 7.2], this implies that M is 
a positive definite matrix and all its eigenvalues are positive. It means that 
eigenvalues of T must be fulfilled A > h. Now, We set A = (T+ hLJy) and 
x =y—Y. Then, using the Lemma 1 for (15), we conclude that 
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||El|2 


Vi. 
lly—- Yl. < min (T + hLJu) 


(19) 
Since T is a nonsingular positive semidefinite tridiagonal matrix and using 
Remark 2, we can easily deduce that omin (T + ALJv) ~ Omin (T). 


u(y(t)) 


Theorem 1. Let the fourth-order derivatives of functions f(t), v(¢, s), 
{ 5 oe at ). 


are continuous. Then for all a € (0,1), we have ||y — Y||, = O(h™™ 


Proof. From the continuity of third and fourth order derivative of the corre- 
sponding functions and according to the relations (9) and (18), we conclude 
that there exists constant ¢ € R such that 


||B||, < ch? + cht*2, 
and so, from (19), we get 


ch? +ch9t2 | th? +chit2 


a 1 
= ch? hi 2, 
Amin (T) > h ee ° 


lly- Yllo< 


According to Theorem 1, the maximum order of convergence is achieved 
when gq > 4. Therefore, we use a Simpson quadrature rule to discretize the 
integral parts of (6). 


Remark 3. As is well known, all norms are equivalent for every z € R”; that 
is, ||Z||1 < Vn|lz|l2 < nlz|oo. As a result, if g > 4, then ||y — Y||, = O(h?) 
and |ly — Y||,, = O(h*) for the proposed method (11). 


3 Numerical examples 


The performance of the proposed method to solve the SFIDE (1) is demon- 
strated in this section. In the following numerical simulations, the step size is 
selected as h = 2-*,k = 2,...,8, and then the error ||y — Y||, is computed. 


Example 1. As a first example for SFIDE (1), we consider u(y) = exp(—y) 
and 


u(t, 8) =u9 t'+9%s?t% sin (unt? *%) sin (urs? tT?) , 
f(@) =t'*7* (3 + a) (—p?0788 4° (3 + a) cos (unt? t*) 
—2 (pum +1)(1+a)sin (unt?**)) 


2Qum(atl1)(a+3)? 


Then the exact solution is y(t) = cos (wt°**), where vp = 
and a € (0,1). 


exp(— cos(uum))—exp(—1) 
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Table 1: Lz error and order of the method (11) for Example 1. 


C= =s Order @=F,p=1 Order a= 7p= 2 Order 
8.9150e-02 — 6.1400e-01 — 3.0755e+00 — 

5.9426e-03 3.90 7.0005e-02 3.18 7.0548¢e-01 2.12 
4.6722c-04 3.66 6.2258e-03 3.49 6.8774e-02 3.35 
3.2741e-05 3.83 4.3233e-04 3.84 8.1641e-03 3.07 
1.1031e-04 3.60 3.5398¢-05 3.61 6.9052e-04 3.56 
128 2.6914e-06 3.52 3.0591e-06 3.53 6.1420e-05 3.49 
256 2.0531e-08 3.50 2.6884e-07 3.50 5.4073e-06 3.50 


RVavrarlse 


It should be mentioned that, in this example, the Jacobian of the kernel 
u(y) is not positive and that increasing 4 causes more oscillations of the 
solution y. We computed the numerical solution of this singular problem 
by the proposed method (11). The numerical results of this test problem 
in the form of the Ly error and the order of the method are reported in 
Table 1. From this table, we can see that the desired order of convergence 
of the presented method is obtained. In Figure l(a), the exact solution of 
the problem given by Example 1 is compared with the approximate solution 
derived by the proposed method (11) when h = 2-8, a = 3, and p = §. 
The absolute error of the present method to solve this example is plotted in 
Figure 1(b). 


-—— Approximate solution 
---- Exact solution 10 


Absolute error 


i— 

-_ uw i—) 
ae 
ie 
i—) 

S 


0.6 0.8 1 


(a) Exact and approximate solutions for h = 2. b) Absolute error. 


Figure 1: Numerical results of the proposed method to solve Example 1 with a = 3, L= 
7 


3: 


Example 2. As a second example for SFIDE (1), we consider u(y) = exp(y) 
and 


u(t, 8) = po t'*t?%s?*° cos (t2*) cos (s*t), f(t) = —(a + 8)? 49 sin (79) . 


Then the exact solution will be y(t) = sin (t°*°), where po = 20s aCe) 


for a € (0,1). 


This example and Example 1 are similar, but the sign of its Jacobian Jy 
of the kernel u(y) is unlike that of the ones in Example 1. We report the 
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1} [Approximate solution 
|---~ Exact solution 10° 
0.8 Fa . as 
5 ce 
a s 
_ 0.6 es 3 Ye \ 
S 2 Neo A 
Sai z 197° i : 
g Me 
2 ut 4 
0.2 < u 
Y 
. 1971 
0 02 04,06 08 1 02 04,06 08 1 
(a) Exact and approximate solutions for h = 278. (b) Absolute error. 


Figure 2: Numerical results of the proposed method to solve Example 2 with a = 3. 


Table 2: Le error and order of the method (11) for Example 2. 


N a= i Order a= 5 Order a=2 Order 
4  2,2839e-02 — 3.1452e-02 — 4.3256e-02 — 

8 2.6238e-03 3.12 3.8926c-03 3.01 5.7418e-03 2.91 
16 
32 


2.1676e-04 3.59 3.3580e-04 3.53 5.1593e-04 3.47 
1.7916e-05 3.59 2.8937e-05 3.53 4.5736e-05 3.49 
64 1.4922e-06 3.58 2.5547e-06 3.50 4.1658e-06 3.45 
128) 1.2372e-07 3.59 2.1784e-07 3.55 3.7527e-07 = 3.47 
256 1.0570e-08 3.54 1.8648e-08 3.54 3.3067e-08 3.50 


numerical results of this example with various values of a and step size h. 
Again, from Table 2, we can observe that the expected order of convergence 
f is achieved. For a = s and h = 2~8, the exact and approximate solutions 
are depicted in Figure 2(a). The absolute error plotting in Figure 2(b) shows 
that the present method is accurate and successful. 


Example 3. As a third example for SFIDE (1), we consider u(y) = —y® and 


1 


u(t,s) =p B2stte (144242)? F(t) = fy 248¢ (1g ete) PPP 


Then the exact solution will be y(t) = (1+ gota)? where vp = suse Oates 
and fo = 6(6 —1)(2+ a)? for a € (0,1) and B > 0. 


We solved this singular boundary value problem for some selected values 
of aand 2. The Lz norm of the errors is computed for the presented method 
(11) and then is exhibited in Table 3. It can be seen that the numerical results 
verify that the claimed order of the convergence of the method is achieved. 


Example 4. Consider the nonlinear SPFIDE from [6] as 


—ey"(t) + (2 — exp(—t))y(t) 4 5 fo (exp(t cos(as)) — 1) y(s)ds = a 
y(0)=1 y(1)=0, 
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Table 3: Lz error and order of the method (11) for Example 3. 


N a=0.25,8=3 Order a=0.5,8=3 Order a=0.75,8=4 Order 
4 9.9948e-01 = 1.5924e+00 = 1.0037e+01 — 
8 4.9934e-01 1.00 2.5881e-01 2.62 2.1306e+00 2.23 
16 7.3623¢e-02 2.76 3.0132e-02 3.10 2.9483e-01 2.85 
32 5.6781e-03 3.69 2.908 1e-03 3.37 3.101 7e-02 3.24 
64 5.0210e-04 3.49 2.6321e-04 3.46 2.8995e-03 3.41 
128 4.0819e-05 3.62 2.3407e-05 3.49 2.635 1e-04 3.45 
256 3.4354e-06 3.57 2.0688e-06 3.50 2.3132e-05 3.50 


where € € (0,1] is a perturbation parameter. Since the exact solution to 
this problem is unknown, we can use the double-mesh principle to estimate 
the errors and compute numerical solutions [6]. The errors obtained so are 
denoted by 

EP = max ro - ye F 
in which 2a is the corresponding approximate solution with respect to ¢ 
and step size h. 


In the reported numerical results, we try to compute the estimated con- 


vergence rates 
Ee EP 
h Pn 
oe (ar) ae (a) : 


where EB” = max, pe, The maximum pointwise errors and the rates of 
convergence p?, De are obtained for the values « = 44,7 = 0,...,4, and 


N = 2',1 = 5,...,10, by the proposed method. In this example, the nu- 
merical results of the presented method will be compared with those of the 
numerical reports in [6]. According to the computational results in Table 
4, we observe that the presented method is more accurate than that of the 
method presented in [6]. It can be seen that the numerical results confirm 
that the methods have achieved the declared order of convergence. Based on 
Table 4, we can conclude that the order of convergence of the present method 
is 4 in the ZL. norm, while the method of [6] is of order 2 with respect to the 
Loo Norm. 


Example 5. As a final example, we consider the following linear singularly 
perturbed boundary value problem [21, 16, 13] 


y(0) = 


< 

— 
= 

Ww 
I 
— 


eo gerbes 0<t<1, 
0, 


The analytical solution of this problem is 
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Table 4: Loo error and order of the methods for Example 4. 


Present method with a= le —6 


E N=2° N = 26 N=2? N=28 N =29 N =210 
1 7.024e-07 4.390e-08 2.751e-09 1.719e-10 1.074e-11  6.528e-13 


3.999 3.997 4.000 4.001 4.040 

4-+  2.591e-06 1.594e-07 9.958e-09 6.218e-10 3.888e-11 2.323e-12 
4.023 4.001 4.001 3.999 4.065 

4-2 3.537e-05  2.638e-06 1.727e-07 1.090e-08 6.827e-10 4.264e-11 
3.745 3.933 3.986 3.997 4.001 

4-3 2.305e-04 1.931e-05 1.283e-06 8.137e-08 5.104e-09 3.197e-10 
3.577 3.912 3.979 3.995 3.999 

4-4 4.185e-04 3.902e-05 2.625e-06 1.669e-07 1.048e-08 6.563e-10 
3.423 3.894 3.975 3.993 3.997 

E” 4,185e-04 3.902e-05 2.625e-06 1.669e-07 1.048e-08 6.563e-10 

pe 3.423 3.894 3.975 3.993 3.997 

Method of [6] 


€ N=2° N =26 N=27 N=28 N =29 N = 210 
1 2.882e-02 7.291e-03 1.839e-03 4.624e-04 1.16le-04 2.906e-05 


1.983 1.987 1.992 1.994 1.998 

4-1 2.86le-02 7.251e-03 1.832¢-03 4.614e-04 1.507e-04 3.785e-05 
1.98 1.985 1.989 1.991 1.993 

4-2 4.001e-02  1.016ce-02 2.575e-03 6.508e-04 1.643e-04 4.139e-05 
1.978 1.98 1.984 1.986 1.989 

4-3 4.331e-02 1.100e-02 2.791e-03 7.070c-04 1.790e-04  4.527e-05 
1.977 1.979 1.981 1.982 1.983 

4-4 4.343e-02 1.105e-02 2.804e-03 7.123e-04 1.809e-04 4.593e-05 
1.975 1.978 1.977 1.977 1.978 

EP  4,343e-02 1.105e-02 2.804e-03 7.123e-04 1.809e-04  4.593e-05 

ph 1.975 1.978 1.977 1.977 1.978 


(2e — 1)(1 — exp(—t/e)) 
1 — exp(—1/e) ; 


y(t) =t(t + 1—2e) + 


In this example, we consider the methods [21, 16, 13] to compare the 
obtained numerical results with the present method. Table 5 contains the 
computed numerical solution achieved by our method and other methods. 
From this table, it can be seen that the present method is more accurate 
than methods [21, 16, 13]. 
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Table 5: Comparison of the approximate solutions of Example 5 for « = h = 1073. 


t Method [21] Method [16] Method [13] Present method Analytical solution 


(a = le — 6) 

0.001 = -1.0009970 -0.6311195 -0.6293169 -0.6298615715 -0.6298573177 
0.010 = -0.9918800 -0.9898546 -0.9878740 -0.9878746043 -0.9878746909 
0.020 = -0.9815600 -0.9796000 -0.9776400 -0.9776399939 -0.9776399980 
0.030 -0.9710400 -0.9691000 -0.9671600 -0.9671599962 -0.9671599999 
0.040 = -0.9603199 -0.9584000 -0.9564800 -0.9564800000 -0.9564800000 
0.050 = -0.9493999 -0.9475000 -0.9456000 -0.9456000079 -0.9456000000 
0.100 = -0.8918000 -0.8900000 -0.8882000 -0.8882000032 -0.8882000000 
0.300 —-0.6114000 -0.6100000 -0.6086000 -0.6086000008 -0.6086000000 
0.500 = -0.2510000 -0.2500000 -0.2490000 -0.2490000008 -0.2490000000 
0.700  0.1894000 0.1900000 0.1906000 0.1906000004 0.1906000000 
0.900 —0.7098000 0.7099999 0.7102000 0.7102000001 0.7102000000 
1.000 = 1.0000000 1.0000000 1.0000000 1.0000000 1.0000000 


4 Conclusions 


In this work, an effective and high-order numerical method for solving a class 
of nonlinear singular two-point boundary value Fredholm integro-differential 
equations was presented. After formulation of the method, as well as uti- 
lizing an appropriate numerical integration, the original problem was con- 
verted to a nonlinear algebraic system. The error analysis was performed to 
demonstrate the robustness of the method. It was observed that the present 
methods achieved the order of convergence O(h™{3-4-2}) in the Lz norm, 
where q is the order of the quadrature method. Here, some test problems of 
type SFIDEs and SPFIDEs are solved numerically. Numerical simulations 
confirmed the theoretical analysis and efficiency of the new method. 
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